Jacobi-like bar mode instability of relativistic rotating bodies

We perform some numerical study of the secular triaxial instability of rigidly rotating homogeneous fluid bodies in general relativity. In the Newtonian limit, this instability arises at the bifurcation point between the Maclaurin and Jacobi sequences. It can be driven in astrophysical systems by viscous dissipation. We locate the onset of instability along several constant baryon mass sequences of uniformly rotating axisymmetric bodies for compaction parameter $M/R = 0-0.275$. We find that general relativity weakens the Jacobi like bar mode instability, but the stabilizing effect is not very strong. According to our analysis the critical value of the ratio of the kinetic energy to the absolute value of the gravitational potential energy $(T/|W|)_{\rm crit}$ for compaction parameter as high as 0.275 is only 30% higher than the Newtonian value. The critical value of the eccentricity depends very weakly on the degree of relativity and for $M/R=0.275$ is only 2% larger than the Newtonian value at the onset for the secular bar mode instability. We compare our numerical results with recent analytical investigations based on the post-Newtonian expansion.Comment: 15 pages, 8 figures, submitted to Phys. Rev.

Two physical characteristics of numerical apparent horizons

This article translates some recent results on quasilocal horizons into the language of $(3+1)$ general relativity so as to make them more useful to numerical relativists. In particular quantities are described which characterize how quickly an apparent horizon is evolving and how close it is to either equilibrium or extremality.Comment: 6 pages, 2 figures, conference proceedings loosely based on talk given at Theory Canada III (Edmonton, Alberta, 2007). V2: Minor changes in response to referees comments to improve clarity and fix typos. One reference adde

Circular orbits of corotating binary black holes: comparison between analytical and numerical results

We compare recent numerical results, obtained within a ``helical Killing vector'' (HKV) approach, on circular orbits of corotating binary black holes to the analytical predictions made by the effective one body (EOB) method (which has been recently extended to the case of spinning bodies). On the scale of the differences between the results obtained by different numerical methods, we find good agreement between numerical data and analytical predictions for several invariant functions describing the dynamical properties of circular orbits. This agreement is robust against the post-Newtonian accuracy used for the analytical estimates, as well as under choices of resummation method for the EOB ``effective potential'', and gets better as one uses a higher post-Newtonian accuracy. These findings open the way to a significant ``merging'' of analytical and numerical methods, i.e. to matching an EOB-based analytical description of the (early and late) inspiral, up to the beginning of the plunge, to a numerical description of the plunge and merger. We illustrate also the ``flexibility'' of the EOB approach, i.e. the possibility of determining some ``best fit'' values for the analytical parameters by comparison with numerical data.Comment: Minor revisions, accepted for publication in Phys. Rev. D, 19 pages, 6 figure

On the Circular Orbit Approximation for Binary Compact Objects In General Relativity

One often-used approximation in the study of binary compact objects (i.e., black holes and neutron stars) in general relativity is the instantaneously circular orbit assumption. This approximation has been used extensively, from the calculation of innermost circular orbits to the construction of initial data for numerical relativity calculations. While this assumption is inconsistent with generic general relativistic astrophysical inspiral phenomena where the dissipative effects of gravitational radiation cause the separation of the compact objects to decrease in time, it is usually argued that the timescale of this dissipation is much longer than the orbital timescale so that the approximation of circular orbits is valid. Here, we quantitatively analyze this approximation using a post-Newtonian approach that includes terms up to order ({Gm/(rc^2)})^{9/2} for non-spinning particles. By calculating the evolution of equal mass black hole / black hole binary systems starting with circular orbit configurations and comparing them to the more astrophysically relevant quasicircular solutions, we show that a minimum initial separation corresponding to at least 6 (3.5) orbits before plunge is required in order to bound the detection event loss rate in gravitational wave detectors to < 5% (20%). In addition, we show that the detection event loss rate is > 95% for a range of initial separations that include all modern calculations of the innermost circular orbit (ICO).Comment: 10 pages, 12 figures, revtex

Various features of quasiequilibrium sequences of binary neutron stars in general relativity

Quasiequilibrium sequences of binary neutron stars are numerically calculated in the framework of the Isenberg-Wilson-Mathews (IWM) approximation of general relativity. The results are presented for both rotation states of synchronized spins and irrotational motion, the latter being considered as the realistic one for binary neutron stars just prior to the merger. We assume a polytropic equation of state and compute several evolutionary sequences of binary systems composed of different-mass stars as well as identical-mass stars with adiabatic indices gamma=2.5, 2.25, 2, and 1.8. From our results, we propose as a conjecture that if the turning point of binding energy (and total angular momentum) locating the innermost stable circular orbit (ISCO) is found in Newtonian gravity for some value of the adiabatic index gamma_0, that of the ADM mass (and total angular momentum) should exist in the IWM approximation of general relativity for the same value of the adiabatic index.Comment: Text improved, some figures changed or deleted, new table, 38 pages, 31 figures, accepted for publication in Phys. Rev.

Relativistic Models for Binary Neutron Stars with Arbitrary Spins

We introduce a new numerical scheme for solving the initial value problem for quasiequilibrium binary neutron stars allowing for arbitrary spins. The coupled Einstein field equations and equations of relativistic hydrodynamics are solved in the Wilson-Mathews conformal thin sandwich formalism. We construct sequences of circular-orbit binaries of varying separation, keeping the rest mass and circulation constant along each sequence. Solutions are presented for configurations obeying an n=1 polytropic equation of state and spinning parallel and antiparallel to the orbital angular momentum. We treat stars with moderate compaction ((m/R) = 0.14) and high compaction ((m/R) = 0.19). For all but the highest circulation sequences, the spins of the neutron stars increase as the binary separation decreases. Our zero-circulation cases approximate irrotational sequences, for which the spin angular frequencies of the stars increases by 13% (11%) of the orbital frequency for (m/R) = 0.14 ((m/R) = 0.19) by the time the innermost circular orbit is reached. In addition to leaving an imprint on the inspiral gravitational waveform, this spin effect is measurable in the electromagnetic signal if one of the stars is a pulsar visible from Earth.Comment: 21 pages, 14 figures. A few explanatory sentences added and some typos corrected. Accepted for publication in Phys. Rev.

Towards a Realistic Neutron Star Binary Inspiral: Initial Data and Multiple Orbit Evolution in Full General Relativity

This paper reports on our effort in modeling realistic astrophysical neutron star binaries in general relativity. We analyze under what conditions the conformally flat quasiequilibrium (CFQE) approach can generate ``astrophysically relevant'' initial data, by developing an analysis that determines the violation of the CFQE approximation in the evolution of the binary described by the full Einstein theory. We show that the CFQE assumptions significantly violate the Einstein field equations for corotating neutron stars at orbital separations nearly double that of the innermost stable circular orbit (ISCO) separation, thus calling into question the astrophysical relevance of the ISCO determined in the CFQE approach. With the need to start numerical simulations at large orbital separation in mind, we push for stable and long term integrations of the full Einstein equations for the binary neutron star system. We demonstrate the stability of our numerical treatment and analyze the stringent requirements on resolution and size of the computational domain for an accurate simulation of the system.Comment: 22 pages, 18 figures, accepted to Phys. Rev.